60 
3 



Of 



Predicting Ranges for Pulsars' Braking Indices 1 

Nadja S. Magalhaes and Thaysa A. Miranda 

nadj asm@gmail . com 

Federal University of Sao Paulo, DCET 
DCET, Rua Sao Nicolau 210, Diadema, SP 09913-030, Brazil 

and 

Carlos Frajuca 

Federal Institute of Education, Science and Technology of Sao Paulo 
R. Pedro Vicente 625, Sao Paulo, SP 01109-010, Brazil 

ABSTRACT 

The theoretical determination of braking indices of pulsars is still an open 
problem. In this paper we report results of a study concerning such determination 
based on a modification of the canonical model, which admits that pulsars are 
rotating magnetic dipoles, and on data from the seven pulsars with known braking 
indices. In order to test the modified model we predict ranges for the braking 
indices of other pulsars. 

Subject headings: pulsars:general; stars: rot at ion 



1. Introduction 

Pulsars are modeled as rapidly rotating, highly magnetized stars composed mainly of 
neutrons. It is observed that their rotation periods are increasing, implying a decay in their 
rotation frequencies. This spin down is quantified by the braking index, n: 

n = (1) 
ft 2 W 
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where f2 is the pulsar's rotational angular velocity and the dot denotes a time derivative. In 
the canonical model the main time- va rying field responsible f or the loss of rotational energy 
in a pulsar is a magnetic dipole field ( lOstriker fc Gunnlll969l ). 



Several theoretical calculations have been tried to explain the o bserv e d values f or the 
br aking index as can be found , for in stance, in lBlandford fc Romanil (119881 ). iMelatod (119971 ) 



or 



Contopoulos fc Spitkovskyj (120061 ). As far as the authors know, however, no theory was 



developed that satisfactorily explains all the values of the braking indices obtained from 
experimental data. In the cases of the pulsars for which this index was measured it lies 
within the range 0.9 — 2.8 (see Tables [1] and [2]) , while the canonical theory predicts n = 3 in 
the approximation of magnetic dipole field dominance. Improvements on this theory have 
been tried since pulsars were discovered, decades ago . It is known that several factors may 



affect the bra king index ([Manchester &: Taylor! 119771 ; iKaspi et al.l Il994j ; iGlendenningi 11996 
Lorimerll2008l ). 



In this paper we analyze a modification of the canonical model aiming at predicting 
braking indices of pulsars based on those that have already been measured. The modification 
consists basically in allowing for a variable braking index in the equation for the braking 
torque. This implies the introduction of a parameter that balances the units in the equation 
that is essentially related to the star's physical characteristics (mass, radius, magnetic field, 
etc.) 

The prediction of braking indices will be tried using two measured parameters: the 
rotation frequency and its first time derivative. The reason for this choice is the possibility 
of predicting ranges of braking indi ces for a significant number of pulsars since many have 
these parameters measured already (lATNF/ CSIROll2012l ) . The determination of exact values 
of braking indices demands more information, like the frequency's second time derivative or 
the angle between the rotation axis and the magnetic axis of the star, but these will not be 
considered in this work. 

In the next section we present a summary of the canonical model and its modified 
version proposed by us. The following sections provide the results using the modified model, 
an analysis concerning them and concluding remarks. 
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2. Models for Pulsars' Spin-Down 



2.1. The Canonical Model and Typical Values 

In the canonical model the law that governs the decay of the rotation is 



with n = 3 (1 Glendenningl 1 1 9 96[ ) and 



K = 



2fi 2 
3c 3 /' 



The moment of inertia of the pulsar, I, is of the order of MR 2 : 

I = XMR 2 . 

The pulsar's magnetic dipole moment is \i: 

fx = R 3 B sin a, 



(2) 
(3) 
(4) 

(5) 



where R is the radius of the pulsar, B is its surface magnetic field and a is the angle of 
inclination of the magnetic axis to the rotational axis. 



The above equations yield 



\B\ 




(6) 



where we used the fact that Q < 0. Typical values for the constants involved in this equation 
will be described as follows. 

The constant A depends on the shape of the pulsar. If it is modeled as a solid ball of 
radius R and mass M then A = 2/5; for a hollow sphere A = 2/3. Depending on the interior 
structure of the star other values may result for A due to the departure of the shape from 
spherical to ellipsoidal, for instance. In the calculations we will use the approximate value 
A ~ 1/2. 

Recent measurem ents indicate the existence of a neutron star with nearly 2 solar masses 
(IDemorest et al.ll2010r). So far this is an iso lated fact, so it still seems reasonable to use the 
canonical value (ILattimer fc Prakashl 12004 ) M = 1.4M Q in the calculations and we will 
assume this number. 

Theoretical values for the star's radius vary from about 0.6 x 10 6 to 1.4 x 10 6 cm 
(ILattimer fc Prakashl 12004 ) so we will adopt the usual value R = 10 6 cm, and we will choose 



sin 2 a = 1 so that lower limits will be implied. 
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Using the above values the first square root in equation (jSJ) yields: 

S = J^fA^ A = 2-3 x l^cm-^B- 8 / 2 . (7) 
V 2 sin a H 

or S = 2.3 xlO 20 Hz 1 / 2 G. This value will be considered the same for all pulsars investigated 
here. 

Returning to equation (jSJ), since in the canonical model all pulsars have the same braking 
index, n — 3, their surface magnetic fields are different from each other as far as there are 
differences in the rotational behavior of the star, given by Q and Q. For example, among all 
pulsars listed in Tables [1] and [3] only the X-ray pulsars present \B\ > 10 14 G while the others 
present magnetic field intensities 10 times smaller or less. These differences have motivated 
the classification of the high \B\ pulsars as magnetars. In this work we will not focus on 
particular differences in \B\ but on general aspects of the pulsar physics. 

Since actual braking indices are not equal to 3, in the next section we investigate a 
modification of this model. 



2.2. The Modified Model 

We now analyze the following model for the pulsar's spin down: 

Ci = -R{t)n naM , (8) 

with n ac tuai being the experimental value of the pulsar's braking index. The main motivation 
in assuming this relation is that it naturally yields n = n actua i since it varies with each pulsar. 
The object K contains a number of physical characteristics of the pulsar like mass, radius, 
intensity of the magnetic field, perhaps internal densities, temperatures, etc. These quantities 



are e xpected to naturally vary with time due to different processes (IBlandford &: Romani 



19881 ). 

Inspired in the basic physics of the canonical model K is given by: 

K = Kb 2 , (9) 

where K is given by equation and b 2 is a positive function of time whose unit in the CGS 
system is ~Hz 3 ~ nactual . This unit (specifically, the difference in the power) suggests that b 2 
carries information related to the correction that must be made in the canonical model for 
the exact braking index be achieved, but at this point the physical content of b is unknown. 
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Using the same procedure as before we find an expression analogous to equation (jfJJ): 



£= \B\b = ±S\l 7 r^—. (10) 

1 >' factual ^ ' 

The object £ contains more physical information about the star than just the intensity of 
its magnetic field due to the presence of b. In this model \B\ cannot be calculated because 
b is unknown. In what follows we will analyze £ instead of b because the right hand side of 
equation ( TTOj) allows estimates using the constants presented previously (through the use of 
equation ([7])) as well as the experimental values of Q, Q and n actua i (or, equivalently, Q). It 
is the presence of n actua i in this calculation that provides more information on the pulsar's 
physics than the canonical model is capable of. 



3. Results 

For the pulsars that have measured values of n we used equation (JTOj) as described 
above and calculated numerical estimates for |£|, which are listed for the different pulsars 
in the second column of Table [2J covering a range from |£| Crab = 1.5 x 10 13 (for Crab) to 
|£| 1734 = 34 x 10 13 (for PSR J1734-3333) in units of Hz( 3 - n -*-'V 2 G. 

In order to predict the values of braking indices we deduced the following relation from 
equations f[7jl - (JIQjl : 

O I _ I O n actual 



which implies 

log 

^actual 



n 



s 2 /e 



logo ■ < u > 

An interesting fact about this relation is that if |£| is estimated then the braking index can 
also be estimated only from the knowledge of the measured values of u and u since S is 
fixed due to the use of typical values of the star's characteristics. 

We found a way to estimate |£| without needing to know \B\ and b, as follows: we will 
assume that pulsars that have values of vq and z>o close to those of the seven pulsars of Table 
[□would have a value of |£| within the range |£| Crab < |£| < |£li734- ^ i s reasonable to assume 
this as far as v , z> and n actua i do not vary much and S is a constant. 

Therefore, given a pulsar with u and z> as mentioned we use equation (TTTj) to estimate 
a range of values for that pulsar's braking index. To this end in that equation we substitute 



-6 - 



l£lcm& ^ or the minimum value of |£| and |^| 1734 for the maximum. We tested this approach 
with the seven pulsars and in Table [2] we present the results. 

We also did this exercise with four radio pulsars and two anomalous X-ray pulsars 
believed to be magnetars. In Table [3] we present the results for these objects which do not 
have known braking indices so far. This approach predicts that the actual values of these 
pulsars' braking indices should lie within the range between n m i nimum and n maxirnurn in the 
respective lines of that table. 



4. Discussion 

For the family of pulsars that we are analyzing the spin frequency and its first time 
derivative are in the range 0.8Hz < vn < 31 Hz and — 3 < z> < —0.01 (xlO -10 s~ 2 ) (see 



Table [2]). According to the catalog lATNF / CSIROl ( 120121 ) there are tens of pulsar within 



these ranges. Using this family the quantity |£| is expected to lie within a range as well, 
given by 1.5 < |f| < 34(xl0 13 Hz( 3 " n -^/ 2 G)(see Table EJ. 

The object £ may vary with time since it depends, for instance, on z/ and z> as well as 
on B (see equatio n (ITuD). For our study , howe ver, we did not need values for the magnetic 



field (a topic that iBlandford fc Romanil (119881 ) explore, for instance). The values of B that 



the pulsars may have will depend on the object b which we did not need to address in 
the present analysis and is open for now. This quantity, like £, is expected to contain 
information on particularities of the star besides its spin behavior, like internal structure, 
stresses, interaction with the environment etc. 

In Table |2] the ranges from u n es timated minimum" to "n estimated maximum" for all those 
pulsars do contain the respective values of their actual braking indices. There are some 
u n es timated r no J ximum" larger than the canonical value (n = 3), indicating that the Crab's 
value for |£| is too small for such pulsars. 

In Table [3] we chose those pulsars as examples for the following reasons: (i) PSR 
B2334+61, PSR J1418-6058 and PSR J1124-5916 have relatively high values for z> thus be- 
ing perhaps good candidates for the experimental determination of i>o; (ii) PSR J1509— 5850 
has a limiting value for z>; and (iii) IE 1547.0—5408 and IE 1841—045 are X-ray pulsars and 
their frequencies or frequency derivatives are somewhat departed from the ranges of interest. 

The first four pulsars in that table show braking index ranges similar to those in Table 
[2j In particular, PSR J1509— 5850 may have a braking index significantly lower than 3. 



The X-ray pulsars show very different values for n es timated maximum which indicate that 



-7- 



they might have important physical differences compared to the pulsars in Table [2j It seems 
that IE 1547.0—5408 has braking index not far from 3. 

5. Concluding remarks 

We analyzed a modification of the canonical model for pulsars' spin down that introduces 
in the equation parameters (£ or b) with physical content of the star. We investigated 
particularly the parameter £ which is sensitive to the star's braking index and thus to its 
physical characteristics. 

Using in the analyzed model the values for the braking index of seven pulsars, obtained 
so far from experimental data, we found ranges for these pulsars' braking indices in which 
all actual braking indices laid. We calculated |£| for these pulsars and used them to predict 
ranges of braking indices for other pulsars, shown in Table [3J These results were possible 
assuming that all pulsars have the same (typical) values for some physical characteristics, 
like mass, radius etc. 

In order to improve the model it is necessary to find physical details about £; this 
would help limit it to better values. It is expected to provide information on the pulsar's 
magnetic field as the definition shows, equation ( TTUl) . but contains also a contribution from 
the unknown object b. We are working on the proposal of a physical context that would 
yield the parameters £ and b and in which mass, radius and other parameters are allowed to 
vary. 

NSM acknowledges Dr. Jose A. S. de Lima for introducing her to the braking index 
problem. TAM acknowledges the Brazilian institution UNIFESP for financial support. NSM 
has a fellowship from the Brazilian federal funding agency CNPq (grant 309295/2009-2). 
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Table 1. Rotation frequencies (z/ ) and their first time derivatives for the pulsars with 
known braking indices. Crab is PSR B0531+21 and Vela is PSR B0833-45. 



Pulsar 




u (Hz) 


z> (xl0- 10 s- 2 ) 


i> (xl0- 21 s~ 3 ) 


Refs 


Crab 




30.22543701 


-3.862283 


12.4265 


1 


PSR B1509 


-58 


6.633598804 


-0.675801754 


1.95671 


2 


PSR B0540 


-69 


19.8344965 


-1.88383 


3.81 


2,3 


PSR J1119- 


6127 


2.4512027814 


-0.2415507 


0.6389 


4 


Vela 




11.2 


- 0.157 


0.031 


5 


PSR J1846- 


0258 


3.0621185502 


-0.6664350 


3.13 


6 


PSR J1734- 


3333 


0.855182765 


-0.0116702 


0.0028 


7 



References. (1) 


Lvne. Pritchard & Graham Smith 


); 


(2) Livingstone et al 


(2007 


): (3) Bovd et al. 


(1995 


Waltevrede. Johnston & EsDinoza 


(2011); (5) 


Lvne et al. 


(1996) 


Livingstone et al. (boil 


): (7) Espinoza et al. (2011) 





(L1993J); 

(4) 
(6) 
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Table 2. Ranges of Braking Indices (n) Obtained with our Approach and the Actual Ones. 



Pulsar 


lei 


^estimated 


^estimated 


n 




xlO 13 


minimum 


maximum 


actual 


Crab 


1.5 


1.3 


2.5 


2.509 


PSR B1509-58 


2.3 


1.4 


3.1 


2.837 


PSR B0540-69 


4.4 


1.3 


2.6 


2.140 


PSR J1119-6127 


7.1 


1.5 


3.8 


2.684 


Vela 


11 


0.9 


2.3 


1.4 


PSR J1846-0258 


19 


1.8 


3.9 


2.16 


PSR J1734-3333 


34 


0.9 


4.6 


0.9 



Note. — In the second column the variable £, detailed in the 
text, is in CGS units. Columns three and four define values for 
the braking indices as estimated by the model; the maximum 
value was calculated adopting the Crab's value for £ (the small- 
est) while the minimum value was obtained using the value of £ 
for PSR J1734— 3333 (the largest). In the column nactual the 
values of the pulsars' braking indices derived from experimental 
data are presented, obtained from the same references shown in 
the preceding table. 
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Table 3. Values for the braking indices of pulsars calculated using the proposed model. 
The maximum value was calculated adopting the Crab's value for £ while the minimum 
value was obtained using the value of £ for PSR J1734— 3333. The sources for the values of 
Uq and z>o are listed in the last column. Pulsar IE 1841—045 is also known as PSR 
J1841-0456 while IE 1547.0-5408 is the same as PSR J1550-5418. 



Pulsar 


u (Hz) 


z> (xl0- 10 s- 2 ) 


^minimum 


^maximum 


Refs 


PSR B2334+61 


2.02 


-0.0079 


0.3 


2.7 


1 


PSR J1509-5850 


11.2 


-0.0116 


0.3 


1.7 


2 


PSR J1418-6058 


9.04 


-0.1385 


0.9 


2.4 


3 


PSR J1124-5916 


7.38 


-0.4100 


1.2 


2.8 


4 


X-ray pulsars: 












IE 1547.0-5408 


0.483 


-0.0541 


2.4 


8.0 


5 


IE 1841-045 


0.848 


-3.216 


1.2 


-9.7 


6,7 



Ref erences. 
f 2009h ; (4) 
H); (7) 



Yuan 



Camilo et al. 



et al. f2nKramer et al.l (l2QQ3h; (3)|Abdo et all 



Kuiper et al. I ( 120061 ). 



2002h: f5) ICamilo et~aD f|2007h : (6) IVasisht fc Gottheli 



